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Abstract 

A rotated and harmonically trapped Bose gas with attractive interactions is expected to either 
remain stationary or escape from the trap. Here we report that, on the contrary, in an anharmonic 
trapping potential the Bose gas with attractive interactions responds to external rotation very 
differently, namely through center-of-mass motion or by formation of vortices, depending on the 
strengths of the interaction and the anharmonic term in the trapping potential. 



*Present address: Department of Physics, KTH, SE-10691 Stockhohn 



1 



Trapped ultracold gases of atoms with attractive s-wave interactions offer quite different 
physics compared to the more common case of repulsive interactions Q]. The behavior under 
rotation is an important example. While repulsive Bose-Einstein condensed gases respond 
to rotation by forming quantized vortices, such states are not necessarily energetically favor- 
able in attracting gases. In fact, Wilkin et al. have studied a harmonically trapped gas of 
attracting bosons under rotation in the limit where the interactions can be treated perturba- 
tively. They find that in the lowest-energy many-body state for a given angular momentum 
(the yrast state) the center of mass (CM) carries all the angular momentum. According to 
Wilkin et al, for this state the single-particle density matrix has more than one eigenvalue 
of macroscopic size, which means that more than one single-particle state is macroscopically 
occupied. Such a state is called a fragmented condensate |3l. The form of the lowest-energy 
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state has since been proved in greater generality Pethick and Pitaevskii [5|] made the 
important observation that the rotational state in question can be seen as a Bose-Einstein 
condensate if it is viewed in the CM reference frame. Nevertheless, it seems that this in- 
teresting state of motion is never thermodynamically stable, because the critical frequency 
is equal to the frequency for destabilization of the cloud. On the contrary, in this Letter 
we show that an anharmonic trapping potential makes the state of CM rotation, as well as 
vortex states, attainable. 

We consider a trapped, attractively interacting Bose gas at zero temperature, subject to a 
rotational force with a given angular velocity, and study the states that minimize the energy 
under such conditions. To this end, we first derive an equation of motion governing such a 
gas, taking into account the CM motion. We then study numerically and variationally the 
behavior in an anharmonic trapping potential, and show that the parameter space is divided 
into regions of CM motion and vortex formation. 

The system under study consists of bosons subject to the many-body Hamiltonian 



The external potential V^(r) can in experiment be of either magnetic or optical origin, and 
is usually taken to be a harmonic-oscillator well. But the harmonic potential has special 
properties that simplify the problem but in fact prevent rotation, and we therefore study 
the more general case of a cylindrically sjnnmetric trap that is anharmonic in the radial 
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direction: 



V{r, 9, z) = -m 
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The oscillator length Oosc is defined as aosc = {h/mtuY^'^. The third term in the Hamiltonian 
is the inter-atomic interaction potential which is here replaced by its usual s-wave ap- 
proximation; the strength Uq = Anfi^a/m depends on the scattering length a which in this 
study will be assumed negative, implying an attractive interaction. 

For the case of purely harmonic trapping (A = 0), the state of CM rotation was discussed 
in Refs. j^, 0, |^ . The energy of L quanta of angular momentum carried by the CM is exactly 
Lhuj, independent of interactions (,7], Appendix A; cf. 0]). This has the important 
consequence that the gas cannot be rotated, because the critical frequency for the excitation 
of any rotation is equal to the frequency where the cloud is destabilized. The energy of a 
system in CM rotation as seen in a frame rotating with angular velocity Q is, namely, 

E -nL, = Eo + n{iu - fi)L, (3) 

where Eq is the energy of the cloud in the absence of rotation. The energy, Eq. (jH]), is 
minimized for L = when Q < uj and for L ^ oo when Q > uj. Consequently, a harmonically 
trapped attracting Bose gas cannot be set into steady rotational motion: it either stays still, 
for slower stirring, or escapes from the trap for faster stirring. 

In order to stabilize a finite-angular momentum state one may add an anharmonic term 
to the trapping potential. In such a potential, the single-particle energy as a function of 
angular momentum increases faster than linear, and the energy in a rotating frame may 
therefore have a minimum for a finite value of L. However, there is a priori no guarantee 
that the lowest-energy state is that of CM rotation when the trap is anharmonic. In fact, 
for a noninteracting gas in an anharmonic trap it is known that the lowest-energy state for 
given angular momentum Lz = Nmih is the state where all particles are in the same angular 
momentum eigenstate^Ji other words, a Bose-Einstein condensate containing a vortex with 
quantum number mi |lO|. This is so because the nonlinear dependence of the energy on 
angular momentum implies that an angular-momentum eigenstate has lower energy than 
any superposition of different states with the same expectation value of angular momentum. 
An intuitive picture can be constructed as follows: in a frame rotating with a fixed angular 
velocity Q, the sum of the centrifugal and potential energies create a Mexican-hat shaped 
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potential if the angular velocity is high enough. In the absence of interactions, the lowest- 
energy state is one of toroidal shape lying along the minimum of the effective potential. When 
the interactions are attractive, however, the tendency for the density to be maximized may 
favor the state of CM rotation. It remains to investigate the crossover between these two 
types of rotational state. 

We now put up an ansatz that allows for the derivation of a simple set of equations of 
motion of an attractive gas. In a harmonically confined gas with isotropic interactions, the 
CM and relative motions decouple and therefore the wave function of the A^-body system 
can be written in the form 

^(ri, . . . , r^) = M'R)M^'i, • • • , r'^), (4) 
where R = (ri + . . . + rM)/N is the CM coordinate and r' = Tj — R are the coordinates 

nn 

of each particle relative to the CM, i = 1, . . . , p, Furthermore, the relative part iIjr 
exhibits Bose-Einstein condensation, so that the associated density matrix has one macro- 
scopic eigenvalue j^. Consequently, the many-body wave function can be written on the 
form 

N 

vl/(ri,...,r^)=^c(R)nx(r:), (5) 

i=l 

as long as the depletion of the condensate is neglected. Although Eqs. (jllE)) hold strictly 
only in the case of a harmonic trap, we shall assume that the decoupling between CM and 
internal motion holds approximately also in the anharmonic case and use Eq. (jSJ as an 
ansatz for the derivation of an equation of motion that can be used to study the rotational 
properties of the attractive gas. We expect the resulting equation to be approximately valid 
when A is small so that the trap does not deviate much from harmonic. 

The interaction energy is independent of the center-of-mass motion, simply because 

^i-T^j = r'i-r'j. (6) 

The harmonic part of the trap potential energy can be expressed in the relative coordinates 
according to 

Ex? = M2 + E^f> (7) 

i i 

and equivalently for the other Cartesian components. The results (jHHZj) together imply that 
in the harmonic case, the CM motion is decoupled from the internal motion, as discussed 
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above. The quartic part of the potential expressed in CM and relative coordinates is 

N N 

Y: rf = NR^ + Y: [2i?Vf + A{Xx[ + Yy[f + 4rf (Xx^ + Yy[) + rf ] . (8) 

1=1 i=l 

The equations of motion are to be obtained in a variational procedure, slightly general- 
izing the derivation given in Refs. ll 11|. We express the energy of the A^-body system in 



terms of the trial wave function, Eq. (0), 

E = jdri... cir^^*(ri, . . . , r^)if^(ri, . . . , r^v), (9) 

where H is the Hamiltonian, Eq. (Q). The transformation to the new system of + 1 
coordinates, R, r'^, . . . , r^, introduces a Jacobian J = N"^ in d dimensions and a delta 
function 6{r[ + r'g + . . . + r^) in the integral, because the relative coordinates are not 
independent. Neglecting this interdependence of the coordinates introduces only errors of 
order A^~^, however, so that the + 1 variables can in effect be regarded as independent 
and the functions tpc and x can be assumed to be individually normalized to unity. The 
trap potential term is rewritten according to Eqs. (0111), and the action of the kinetic-energy 
operator on the trial wave function is readily calculated. With a number of superfluous 
factors x(r9 integrated out, the energy functional is 

92 Nm 



E = J JdRdr''^*c(R)x*{r' 
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(r'' + 2R\'^ + A{Xx' + Yy'f + Ar'\Xx' + Yy' 



^c(R)x(r'), (10) 



where J = as noted above. The energy is now functionally minimized with respect to 
the two wave functions ipc and x- The result is the two coupled equations 
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+ Vc{R) ^c(R) = EcM'R), (11) 



Vn{r') + Uomr)\' } $(r) = /i$(r). (12) 



[ 2m dr"^ 

We have introduced here the condensate wave function $ = ^/N x and dropped the prime on 
the relative coordinate r. The effective potentials for the center-of-mass and relative motion 
are defined as 

VciR) = Imu;\1 + 4A^)i?2 + ^R\ (13) 



Vn{r) = l-muj^l + AX^-^y + (14) 

^ '^osc ^'^osc 

where (i?^) = / JdIlR'^\ipci'R)\'^ and (r^) = N^^ J (irr^|$(r)p. Averages of odd powers in 
the relative coordinate vanish because the relative coordinates are centered around the CM. 
The quantities Ec and /i are introduced as Lagrange multipliers to ensure separate normal- 
ization; Ec is actually the center-of-mass energy and /i is the chemical potential as usual, 
and M = Nm is the total mass of the gas. For harmonic trapping (A = 0), the equations de- 
couple: the CM behaves like a single particle of mass M, and the condensate wave function 
for the internal motion obeys the Gross-Pitaevskii equation. The same effective potential 
for the internal motion was derived in Ref. 3]; however, the Gross-Pitaevskii equation was 
there taken starting point. 

In the absence of rotation, the spatial extent of the CM wave function is on the order 
of the effective CM oscillator length ac = {h/McuY^'^ while the size of the condensate is 
of order Oosc, which is a factor A^^^^ larger. As a result, the CM motion has a negligible 
effect on the motion of the condensate, and the usual Gross-Pitaevskii equation is retained. 
However, we are interested in the case of highly excited CM motion, where the total angular 
momentum L carried by the CM is comparable to the number of particles A^, and in that 
limit, the effects of the CM wave function cannot be neglected. 

Numerical solutions to the coupled equations of motion ()1HI14|1 are obtained with the 
steepest descent method. The z dependence of the CM wave function ipc decouples and 
is exactly Gaussian, while for the condensate wave function $ both the radial and axial 
dependence is subject to numerical minimization. We have chosen the harmonic part of the 
potential to be isotropic, Uz = oj. For large values of L, the CM wave function becomes very 
peaked and difficult to handle numerically. It has therefore been necessary to restrict the 
number of particles N to N < 15, and the possible angular momenta L to a few times A^. 

The result for the total angular momentum is shown in Fig.^ We have chosen Na/aosc = 
—0.1, and A = 0.05 and A = 0.1 as two illustrative examples. Both these parameter values 
lie in the region where for slow rotation the CM carries the motion, hence the total angular 
momentum per particle L^/N increases only in steps of hN~^, but for faster rotation the 
angular momentum is carried by a vortex in the condensate wave function, and the angular 
momentum per particle is integer. For sufficiently large A and weak attraction, vortex states 
are energetically favorable for all Q, because the anharmonic potential favors vortex states 
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where the density distribution is toroidaL Conversely, for large enough Na/aosc and small 
enough A, the attraction between bosons favors the state of CM rotation where the density 
distribution is maximally concentrated. 

The numerical results give at hand that a mixed state, where both m/ and L are nonzero, 
is never the energetically favorable situation. This is clear, since the density distribution 
of such a state is a toroid that is not centered at the middle of the trap, but occupies 
regions of high potential energy in the Mexican-hat shaped effective potential formed by the 
anharmonic and centrifugal energies. 

In order to be able to cover larger parameter spaces, we must consider a variational 
solution of the equations of motion (|imi4|) . We choose trial functions with Gaussian envelope 
for the CM and condensate wave functions: 

$(r) = Ar^(a; + ii/)™'e-'''/2"'*°-"("^/^)"'/2('^"°=^)'. (15) 

That this is a good ansatz has been confirmed by inspecting the shape of the numerically 
calculated wave functions. Mc and Mr are normalization constants. As already noted, the 
z dependence of the CM wave function is not subject to variational minimization, but is 



a Gaussian with the fixed width ac^^j^z- On the other hand, the widths 6, c and d are 
variational parameters. The ansatz allows for the angular momentum to be shared between 
the CM and internal motion and thus incorporates both the known limits of a BEG in a 
vortex state for f/g = 0, and CM rotation for A = 0. The total angular momentum is 
= (L + Nmi)h; we have already seen that the lowest-energy states in a rotating frame 
have either L = or m/ = 0. Defining q = L/N, the energy is 

N ^ /I 



E = —qhuj + ^ + + l)A6c + Ag6^^ 



N , 



- + c + X{mi + 2)c^ 
c 
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+Nnu{ , ^ -. (16 

\2tiuJz) aosc (m;!)222™< erf ^ > 

Terms of order unity are discarded. The energy is now minimized with respect to the 
parameters 6, c and rf, and the energetically optimal way of distributing the total angular 
momentum between the CM and internal motion is determined. The result is plotted as 
solid lines in Fig. ^ It should be noted that the variational solution is carried out for the 
limit of large N while the full numerical solution is restricted to finite values of N . 



Fig. 121 shows the boundary between the CM-rotation and vortex regimes for a few fixed 
angular momenta as calculated in the variational scheme. As expected, a large A, large 
angular momentum and weak coupling favors the vortex state. The coupling at which 
the cloud becomes unstable against collapse is also shown; the apparent shrinking of the 
stable regime with increasing A is due to the decrease of the effective size of the cloud, and 
therefore an increase in density, when A is increased at fixed uo. Included is also the result of 
a perturbative calculation, treating the exact many-body wavefunctions to first order in A 
and Na/aosc (cf- 2101); it is seen that perturbation theory only works when these quantities 
are very small. 

In conclusion, a Bose gas with negative 5-wave scattering length confined in an anhar- 
monic trap will responds to a rotational force by either exciting center-of-mass rotational 
motion or forming a multiply quantized vortex, depending on the strength of the interaction 
and the anharmonic term. The wave function of the CM position and the condensate wave 
function associated with the internal motion are described by two equations that in the gen- 
eral anharmonic case are coupled. These results show that the behavior of condensates with 
attractive interactions is markedly dependent on the harmonicity of the potential, and that 
an anharmonic trapping potential is expected to allow for vortex creation. Reference 3| 
reported experiments using an anharmonic potential with a small quartic term (A ~ 0.001), 
created by imposing a Gaussian optical potential over a magnetic trap. A similar setup, 
possibly with a stronger anharmonicity through increased laser power or decreased beam 
width, could be employed in order to explore the rotational states discussed in this paper. 

While preparing this manuscript we became aware of the related work by Ghosh |l4| . 

Discussions with C. J. Pethick are gratefully acknowledged. The authors acknowledge 
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Na/a =-0.1,X=0.1 
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FIG. 1: Dependence of angular momentum on rotational frequency. The coupling Na/uosc and 
anharmonicity A are as indicated above each panel. The thick solid lines indicate the variational 
result, and the thinner lines represent the fully numerical results for the cases A'^ = 15 (full), 10 
(dashed) and 8 (dotted). When the angular momentum is carried by the center of mass, the curves 
display a smooth increase, but when vortices appear in the condensate, the angular momentum is 
fixed to integer values. 
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FIG. 2: The solid lines represent the boundaries between the three phases; vortex phase, center-of- 
mass rotation and unstable phase; for the case of one quantum of angular momentum per particle, 
Lz = Nfi. The dashed lines represent the case = 2Nh. The dotted line is the boundary between 
the stable and unstable regimes for a nonrotating cloud. The lines are calculated in the variational 
scheme. The dot-dashed line is the perturbative result. 
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